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We investigate the excitation of quantum plasmonic states of light in graphene using end-fire and 
prism coupling. In order to model the excitation process quantum mechanically we quantize the 
transverse-electric and transverse-magnetic surface plasmon polariton (SPP) modes in graphene. A 
selection of regimes are then studied that enable the excitation of SPPs by photons and we show 
that efficient coupling of photons to graphene SPPs is possible at the quantum level. Futhermore, 
we study the excitation of quantum states and their propagation under the effects of loss induced 
from the electronic degrees of freedom in the graphene. Here, we investigate whether it is possible 
to protect quantum information using quantum error correction techniques. We find that these 
techniques provide a robust-to-loss method for transferring quantum states of light in graphene over 
large distances. 


I. INTRODUCTION 


Quantum plasmonics is attracting considerable inter¬ 
est at present from a wide range of researchers, most no¬ 
tably those in the quantum optics and plasmonics com¬ 
munities [1]. This interest is in part due to the poten¬ 
tial of plasmonics for applications in quantum informa¬ 
tion processing (QIP), which include ultracompact and 
versatile single-photon sources ms], and single-photon 
switches [6H8]. An important aspect in the quantum 
study of plasmonic systems is the excitation of quan¬ 
tum states of propagating surface plasmon polaritons 
(SPPs) using photons, which are easier to generate ex¬ 
perimentally in a well controlled manner [9]. Indeed, 
the ability to transfer a quantum state between two dif¬ 
ferent systems is an important requirement for QIP in 
general m- Here, the quantum state transfer must 
be efficient and not entail signihcant decoherence. Fur¬ 
thermore, the transferred state must be maintained over 
times and distances necessary to perform the required 
processing operations. Experimental work in quantum 
plasmonics using photons to excite SPPs has so far con¬ 
firmed the preservation of entanglement when transfer¬ 
ring quantum states between photons and SPPs mm, 
the preservation of superposition states [13], as well as 
a wide range of other properties related to the quantum 
statistics of the excitation process [IlHIHI- Most recently 
experiments have demonstrated two-plasmon quantum 
interference, confirming the maintenance of the bosonic 
nature of the photons used to excite SPPs in a plasmonic 
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FIG. 1: Excitation of quantum plasmonic states in graphene, 
(a): Free-standing graphene and coordinate system. Here, 
SPPs cannot be excited by photons directly without a cou¬ 
pling method to provide energy-momentum matching condi¬ 
tions. (b): End-fire method for coupling, (c): Prism method 
for coupling, (d): Using either end-fire or prism methods, the 
excited plasmonic state propagates along the graphene sur¬ 
face. The hybrid nature of the SPP - consisting of a joint 
state of a photon and a collection of electrons - means that 
the electronic degrees of freedom induce loss effects in the 
photonic part. An error-correction code is introduced to deal 
with the loss and provide propagation over a large distance. 


Most work on quantum plasmonic systems has so far 
focused on basic metallic material as the support media 
for the plasmonic excitations, using either silver or gold. 
However, there is a large range of other materials avail¬ 
able to use [23|, and graphene has recently emerged as 
a powerful alternative due to the possibility of chemi- 
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cal doping and electrical gating, which provides a highly 
tuneable media for supporting quantum plasmonic sys¬ 
tems. Here, studies have focused on emitter coupling and 
decay into SPPs [21] , active control over a quantum state 
biasing | 2 S], nonlinear quantum optics [2Sl428j . quantum 
networks [29] and quantum sensing applications [30]. De¬ 
spite some impressive work in this area, one key issue that 
has not been looked at in detail is the efficiency of the 
excitation of single SPPs in graphene using photons and 
more general quantum states of light. Moreover, given 
the many advantages that graphene offers compared to 
conventional plasmonic media in terms of tuneability, the 
excited SPPs still suffer from the effects of loss as they 
propagate. 

In this work we investigate the two issues of the excita¬ 
tion efficiency and the effects of loss by studying the exci¬ 
tation of quantum plasmonic states of light in graphene 
using end-fire and prism coupling methods. We quan¬ 
tize the transverse-electric and transverse-magnetic SPP 
modes in graphene in order to build a fully quantum me¬ 
chanical model for the excitation process. We then study 
various parameter regimes that enable the excitation of 
single SPPs by photons and find that efficient coupling of 
photons to graphene SPPs is possible. Furthermore, we 
study the subsequent propagation of the excited quan¬ 
tum states under the effects of loss induced from the 
electronic degrees of freedom in the graphene. In order 
to protect the quantum states from loss we use a quan¬ 
tum error-correction code and find that the code pro¬ 
vides a robust-to-loss mechanism for propagating quan¬ 
tum states of light in graphene over large distances. 

The work is divided into five sections. In Section II we 
introduce the model for SPP quantization in graphene. 
In Section III we then use this model to study the con¬ 
version of single photons to single graphene SPPs, and in 
Section IV we provide a model for describing the effects 
of loss during the subsequent propagation of the SPPs. 
In Section V we then investigate quantum state transfer 
and propagation in detail, introducing an error correction 
code for protecting against the effects of loss. We show 
the benefits of using the code compared to not using it. 
In Section VI we conclude with a summary of our results 
and an outlook on future studies. 


II. GRAPHENE SPP QUANTIZATION 


In our study of the transfer of quantum states between 
photons and graphene SPPs we consider the graphene as 
a free-standing sheet, as shown in Fig.[^(a). Here, SPPs 
are excited by photons by using either the end of a fibre 
(end-fire method), as shown in Fig. (b), or a prism, 
as shown in Fig. (c). In order to model the coupling 
between photons and SPPs we must first quantize the 
SPPs in graphene. In this section we briefly summarize 
the quantization steps. 


1. Graphene Conductivity 

We start by considering a laterally-infinite graphene 
sheet lying in the x — y plane, as shown in Fig. (a). 
The graphene is modeled as an infinitesimally-thin, local, 
two-sided surface characterized by a surface conductivity 

CT [SI], 


cr(w) 


ie'^ksT I 
'nffl (w -|- iF) \ 
+ ij) 
ttIt? 


Me 

keT 



fd{-e) - fd{e) 

(cj + n)2 - 4(e/fi)2 


( 1 ) 


where w is the radian frequency. Me is the chemical po¬ 
tential (Fermi energy) and F ( 7 ) is a phenomenological 
intraband (interband) scattering rate. In addition, T is 
the temperature, e is the charge of an electron, is 
Boltzmann’s constant and fd(£) = ( 6 ^"^“^=^/^®^ -I-1)“^ is 
the Fermi-Dirac distribution. The first term in Eq. 0 is 
due to intraband contributions, while the second term is 
due to interband contributions. Absorption is associated 
with both intraband electron scattering and interband 
electron transitions. While Eq. Q is valid for arbitrary 
T, when <C (|Mc|, ^), *-e- the low temperature limit, 
it becomes [24l 
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where Q{x) is the Heaviside function. This form of the 
conductivity has no T dependence and is used in our 
study to simplify calculations in the low temperature 
limit, whereas the full form given in Eq. 0 is used 
in the high temperature limit, i.e. room temperature 
(T = 300K). 

For the decay rates, typical intraband scattering times 
are t = 1/F = 0.35 ps at room temperature, and as large 
as T = 3 — 5 ps at low temperature [33H36] . For the inter¬ 
band scattering rate we use 1/7 = 0.0658 ps [37]. These 
values are assumed throughout our work unless otherwise 
noted. The Drude form of the conductivity (the first 
term) has been verified in the far-infrared |34II37j . and 
in the near infrared and visible the interband behavior 
has been verified [3S]. Note that for the present applica¬ 
tion one could also consider N closely-spaced graphene 
monolayers, modeled as a single layer with a larger effec¬ 
tive conductivity, creff(w) = Na{u!). 


8. Classical Vector Potential 

We start by working in the Lorentz gauge (E = 
-dA/dt, B = V X A) and consider a homogeneous ma¬ 
terial having relative permittivity, e^, on either side of 
the graphene, the wave equation for the vector potential 
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is 

A± (r, i) - ^ ^A± (r, t) = 0, (3) 

where ± corresponds to the region and c = Co / 
with Co being the speed of light in vacuum. The structure 
is invariant in the transverse plane and therefore we can 
use the solution 

A±(r,t) = ^ , 

k 

(4) 

where the wave vector k = fc^x + kyy is parallel to the 
interface. Substituting Eq. Q into Eq. (|^ we find 

(^-9o)AJ(^) = 0, (5) 

where Qq = with k^ = |kp. Thus, a solution 

for the vector potential is 

A± (r, t)=Y^ ^ ^ ^ (g) 

k 

Assuming for the moment propagation in the x direction 
only, i.e. invariance in the y-direction (ky = 0), for TE 
modes the non-zero field components are (Ej,, H^) 

arising from Ay, and for TM modes we have non-zero 
components {Ex^ E^, Hy) arising from and A^. En¬ 
forcing the boundary conditions [38] 

z X (H+-H") = J = ctE (7) 

z X (E+ -E") = 0 

at z = 0 leads to the determination of the form of the 
AJ vectors and to the dispersion relations for TE and 
TM graphene modes that link the wavenumber k to the 
frequency ujk via the conductivity, cr{u!k), and the relative 
permittivity of the background material, Sr- The explicit 
form of the A^ vectors is given later in this section. First, 
we discuss the dispersion relations, which are given by 



where ko = Wfe/co, rj = y/ /io/ir/^o^r (with = 1) and 
cr = (T(a;fc). 

By decomposing the conductivity into real and imagi¬ 
nary parts, a = a' ia", it can be shown that that for 
cr" > 0 (inductive surface reactance) only a single TM 
surface plasma wave is supported by the graphene and 
can propagate; if a" < 0 (capacitive surface reactance) 
the TM SPP is on the improper Riemann sheet, expo¬ 
nentially increasing as |z| —>■ oo EillD]. 



Frequency (THz) 



FIG. 2: Graphene conductivity and surface plasma wavenum¬ 
ber. (a): The conductivity (y{ui) for bare graphene. The 
conductivity has been rescaled by (Jmin = 7re^/2/i. (b): The 
rescaled surface plasma wavenumbers for TE and TM modes 
(sr = 1), AC = kx/k(, = k! -\- in". In both, the chemical poten¬ 
tial is /Xc = 0.5 eV, and we have set r = 5psatr = 0K and 
T = 0.35 ps at T = 300 K. 


On the other hand, for cr" < 0 (capacitive surface re¬ 
actance) only a single TE surface plasma wave can prop¬ 
agate (if a" > 0 the TE SPP is on the improper Rie¬ 
mann sheet) (311 HO]. Fig. (a) shows the conductiv¬ 
ity cr and Fig. (b) shows the rescaled surface plasma 
wavenumber, n = k/ko for both TE and TM modes, de¬ 
composed into real part ac' = Re(fca;/fco) and imaginary 
part k" = liJi{kx/kQ) over a wide range of frequencies for 
T = 0 K and T = 300 K, with the chemical potential 
chosen as /Xc = 0.5 eV. At low frequencies huj 2/Xc 
(w = Wfc is used for concise notation) the intraband con¬ 
ductivity is dominant (cr" ~ o’intra)! interband absorption 
is blocked, and a slow TM surface plasma wave can prop¬ 
agate on the graphene surface (cr" >0). As the frequency 
increases the Drude term falls off, and in the vicinity of 
fiu;/ 2 /Xc = 1 interband absorption becomes important. 
As the frequency increases further cr" ~ cr;"^^^., so that 
a loosely-bound TE surface plasma wave can propagate 
(cr" < 0). Therefore, in Fig. 2(b), the mode is TM to 
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the left of the discontinuity (TE mode is on the improper 
Riemann sheet), and TE to the right of the discontinuity 
(TM mode is on the improper Riemann sheet). Note that 
the frequency range of the two modes can be adjusted by 
changing the chemical potential, which will shift the dis¬ 
continuity to the left or right as required. 

The vector potential from Eq. <§ can be written ex¬ 
plicitly as 

A± (r, t)=Y, + C.C., (10) 

k 

where the Ck are mode amplitudes and the mode func¬ 
tions are 


For TE modes (cr" < 0) the energy stored in the 
graphene can be obtained by temporarily assuming the 
graphene has a small but non-zero thickness d with ef¬ 
fective permittivity 




with associated energy 
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Note that, from the relation and 

Eqs. Q and (|^, we have 


where Ay is the component of potential parallel to the 
graphene sheet. Note that at this stage we are not in¬ 
cluding loss, so that cr = ia". 

For TM modes cr" > 0, leading to a negative permit¬ 
tivity e < 0, and so a different method must be used. We 
again assume that the graphene has a small but non-zero 
thickness d. Inside the graphene the equation of motion 
F = msL leads to 
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with the sign chosen so that the modal field decays expo¬ 
nentially in the vertical direction away from the graphene 
sheet. As cr(w) is a complex valued function, in general 
qo = q'o + iq'o, where gp describes the strength of lat¬ 
eral field confinement (the decay length) to the graphene 
sheet while corresponds to the wavelength of the prop¬ 
agating mode in the z-axis, corresponding to leakage into 
the far-field. For simplicity, in this work we focus on 
propagating modes in the x — y plane only, so we re¬ 
strict our interest to the regime where Re(go) Im(go). 
Furthermore, our model will only include loss in a; — y 
plane, corresponding to the case where internal loss is 
much greater than leaky loss originating from gp. 


3. Quantization 


With the explicit form of the vector potential now 
given, we proceed to quantize the surface plasma wave. 
Note that here we are considering an ideal case with no 
damping effects due to the electronic degrees of freedom 
in the graphene sheet. Furthermore, and consistent with 
no damping, for quantization we assume a non-dispersive 
model since we consider quantum states that are associ¬ 
ated with wavepackets that have very narrow bandwidths 
centered on w. Both assumptions simplify the quantiza¬ 
tion procedure. Damping will be reintroduced to the 
model later in Section IV. 

We start with the Hamiltonian for the field given by 



dA 

Ih 


-kcg|V X Ap dV. 


where v is the velocity of electrons in the elec¬ 
tron gas. Thus, A|[ (r,t) = (m/e)v(r,t) with as¬ 
sociated current density J = —en (r, t) v (r, t) = 

— (e^/m) n (r, t) Ay (r, t), where n is the number density. 
The energy stored in the graphene electron gas kinetics 
is [H] 


He = J n (r, t) (r, t) dV 

= / n (r, t) A| (r, t) dV, (18) 

2 m Jy " 

Considering that for a lossless plasma rie = 
WpEo {mje^) = moja"/e^d, we obtain 

H™ = \j^5{z)ura"A\dV. (19) 

Therefore, the total Hamiltonian for either TE or TM 
modes is 
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dV\ e 



( 20 ) 


-kcg|V X A|Vd(z) 


W cr 
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To compute the Hamiltonian, we take a region of 
space of size L x L on the surface of the graphene sheet 
(in the x — y plane) such that the wavenumbers are 
kx,y = 2,'Knx^y/L, rix^y = 0,±1,±2,..., i.e. periodic 
boundary conditions. The volume integrals in the Hamil¬ 
tonian are evaluated using fj^xdy = L^dk.p 
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and substituting in Eq. (10). By doing this we obtain 
terms such as 


\Air,t)\^ = 2j2CkCMr)-Kir) 

k 

+ ^ C'kC_k</)k(r) • 0-k(r)e-‘2-‘‘ 


+ c.c. 


It is straightforward to show that the cross terms, 
e.g. 0k(r) • (/>_k(r), associated with electric, magnetic, 
and gas kinetic energies will cancel, as occurs in free- 
space optics, and thus we only retain terms such as 
</>k(r) • = l</’k(r)P- By collecting all the non-zero 

terms, after substituting in Eq. (10) and carrying out the 
integrals, we obtain 


H = eoL^ Y. + CkCi:), 


( 21 ) 


where = ^te/tm ^ normalization parameter hav¬ 
ing units of length. To quantize the fields, we use the 


correspondence of the Hamiltonian in Eq. (21) with that 


of the harmonic oscillator, mapping the coefficients as 


^ vSS;'”" ^ vSS;'’"' 

where the bosonic creation and annihilation opera¬ 
tors satisfy the commutation relation [fok 7 ^k/] = i^kk'- 
We then have the quantized Hamiltonian, H = 
'^y{hujy/‘i){J}yb\. + &kSk), for the system and the cor¬ 
responding quantized vector potential operator is 


A(r,t) = Y 


h 






z(k-r— 


H.c. 


. 

For our purposes it is convenient to take the continuum 
limit using —J- (^)^ / dk and by —>■ (^)6(k), so that 


. 

For simplicity, we also assume excitations propagating in 
the x-direction only with a beamwidth W in the y-plane. 
Then, f dk^ and b{k) 

which leads to 


H.C., (25) 

with the spatial mode functions <p{z, k^) given in 
Eqs. ([^-(@. Finally, converting to the frequency do¬ 
main using kx = v~^uj, dkx = v~^dio and b(kx) —>■ 
ydi^6(a;), where Vg(u}) = dio/dk^ is the group velocity, 
we have the vector potential for quantized surface plasma 


a) 



Frequency (THz) 


FIG. 3: End-fire coupling of photons to SPPs on graphene 
using a fiber. Inset shows the configuration considered, 
(a): Transmission coefficient /3- (b): Power distribution for 
TE excitation of graphene strip showing good coupling to a 
graphene TE SPP. (c): Power distribution for TM excitation 
of graphene strip showing no excitation of the graphene TM 
SPP. The cladding has radius 3 pm and Sdadding = 1.16, the 
core has radius 500 nm and Score = 1.45, the depth of the 
removed-region is 3.025 pm, and the graphene is biased at 
Pc = 0.02 eV, with T = 0 K, and r = 0.1 ps (Im(CT) < 0, 
allowing only TE mode propagation). 


waves, now surface plasmon polaritons (SPPs), given by 

AsPP ^2£oWUgWAfTE7TM 

X ^™/™( 0 , w)e-*‘^(*-^/''»)6(a;) -b H.c. 


III. PHOTON-TO-SPP COUPLING MODEL 

With the graphene SPPs quantized we now introduce 
the coupling model between photons and SPPs in order 
to investigate the efficiency of the excitation process at 
the quantum level. As described in Refs. [42l|43], within 
a linear response regime the coupling of photons to SPPs 
can be described in the Heisenberg picture by a unitary 
transformation matrix 


dout (^) 

dout (^) 

where | 7 (w)P -I- |/3(w)P = 1 and d{uj) is an annihi¬ 
lation operator for the photon field which, together 
with d^(w), satisfies the bosonic commutation relation 
[d(w), (^^(a;')] = 6{u} — to'). In the following, we will use 
the coupling coefficient, g{oj) = sin“^ I/3(‘^)|7 de¬ 

fined in terms of the transmission coefficient, Piui), from 
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I3{uj) 

di„(w) 
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Eq. (26). The coupling coefficient appears in the in¬ 
teraction Hamiltonian for the system, given by Hint = 
ih J du}[g{uj)a)b — g*{u})ab^. Here, perfect coupling cor¬ 
responds to g{u!) = 7r/2, which provides the complete 
transfer of a given quantum state of a photon to a SPP. 


In appendix I we provide details of various different 
coupling scenarios that can provide a range of transfor¬ 
mation matrices for the photon-to-SPP transfer. Here 
we briefly summarize the main results. The first cou¬ 
pling scenario we consider is end-fire coupling, where the 
photon field from the end of a fiber is evanescently cou¬ 
pled to the field of the graphene SPP . We consider a 
structure similar to that shown in Fig.^(b), where an op¬ 
tical fiber has a center perturbed region with the cladding 
removed. In order to make the setting more realistic we 
consider the core of the fiber as being partially removed in 
order to support the graphene sheet, as shown in the inset 
of Fig. 1^ (a). Here, photons enter the fiber from the far- 
right and couple to the graphene SPPs, after which the 
output is taken on the left hand plane where the structure 
terminates. We use a numerical FDTD simulation (see 
Appendix I for details) to obtain the transmission coeffi¬ 
cient /3, which is shown in Fig.[^(a). The graphene SPPs 



d (nm) 



in this geometry are quantized using the formalism given 
in the previous section, with appropriate consideration 
of the asymmetric dielectric media above (air) and below 
(fiber core support). The FDTD simulation enables the 
calculation of the overlap between the evanescent photon 
modefunction and the graphene SPP modefunction. As 
the modefunctions correspond to the classical wavelike 
part of the photons and the SPPs, the overlap calculation 
is essentially a classical calculation with the resulting (3 
value entering into the transformation matrix of Eq. (26) 
in order to model the coupling quantum mechanically. 
Here, the operators are associated with a given spatial 
modefunction at a specific frequency uj. 

In Fig.[^(a) one can see that at high frequencies (where 
only TE SPPs are supported), good TE SPP coupling 
can be achieved, whereas TM SPP coupling is signifi¬ 
cantly reduced. In Fig. |^(b) and (c) we show the power 
distribution for the transfer of the field from the fiber to 
the graphene sheet for TE and TM modes respectively. 
One can see that the TE mode can be coupled to well, 
whereas the TM mode cannot as it is only supported at 
lower frequencies. In Appendix I we discuss a grating 
method for efficient end-fire coupling of photons to TM 
modes at lower frequencies, which is tunable by varying 
the chemical potential. 

The second coupling scenario we consider is prism 
coupling, where the photon field from below a prism 
is evanescently coupled to the field of the graphene 
SPP PSHTT] . The coupling structure is similar to that 
shown in Fig.[P(c), whose configuration is shown in more 
detail in Fig. [^(a). In Appendix I, we provide details of 
the prism coupling, which, unlike end-fire coupling, we 
are able to obtain an analytical solution easily and there¬ 
fore do not need FDTD simulation. In Fig. (b), as an 
example, we show the coupling coefficient /? as the spac¬ 
ing between the prism and the graphene, d, changes for 
both TE and TM SPPs at a specific frequency for the in¬ 
coming photon, corresponding to a free-space wavelength 
of Aq = 500 nm. In the inset we show the TE reflection 
coefficient near the resonance of the coupling. It can be 
seen from Fig. (b) that good photon-SPP coupling can 
be achieved using a prism for both TM and TE SPPs. 

In summary, both end-fire and prism coupling methods 
can provide good photon-to-SPP coupling, with trans¬ 
mission coefficients /3 > 0.7 {g > 0.77) for TM and TE 
modes. 


FIG. 4: Prism coupling of photons to SPPs on graphene, 
(a): The coupling system considered, (b): Transmission co¬ 
efficient for SPPs as a function of prism-graphene spacing d. 
TE mode (dotted, prism-graphene spacing given by top hor¬ 
izontal axis): 6i = 54.74°, ei = 1.5, g,c = fujj/2 = hc/Xo 
i/ic = 1.24 eV), / = 600 THz (Ao = 0.5 pm), T = 0 K, 
TM mode (dashed, pc = 0.5 eV, solid, pc = 0.8 eV, prism- 
graphene spacing given by bottom horizontal axis): / = 0.81 
THz (Ao = 3.7 X lO"'^ m), 9i = 64°, ei = 1.5, T = 300A. 
The inset shows the TE reflection coefficient as a function of 
frequency near the resonance frequency. 


IV. PROPAGATION AND DAMPING MODEL 

Once the SPP is excited using one of the above meth¬ 
ods it propagates along the graphene surface. Here, it 
is damped by interactions with phonons, impurities, and 
defects at both the light level (diffraction and radiation 
at a physical discontinuity) and electron level (intra¬ 
band electron scattering and interband absorption), as 
well as with interactions with the thermal bath of field 
modes. The former (possible diffraction and radiative 
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scattering) is ignored here as we assume an unperturbed 
graphene surface. On the other hand, both electron- 
level damping (incorporated in the graphene conductiv¬ 
ity in Eq. 0 ) and thermal interactions cannot be ne¬ 
glected and are accommodated by using a standard mul¬ 
tiple beam splitter model [42l |43l IdSUbC)] . consisting of 
N quantum beam splitters each with a quantized field 
mode Ci{uj), i = 1, These bath field operators sat¬ 

isfy the bosonic commutation relations [ci(a;), c|(w')] = 
SijS{uj — oj'). In the continuum limit N ^ oo, Ax —>■ 0, 
Ci(a;) —> VAa;c(a;, x'), and 6ij —)■ Ax5{x—x'), and the an¬ 
nihilation operator of the SPP after travelling a distance 
X is given by [HJ 03] 


5out(a;,x) = e*'=“="Sout(w) (27) 

+ i^y2koK''{u;) [ dx'^'^-^^-^''>c{uj,x'), 

Jo 


where = k^iui) = k^i^n'{ lo) in"( lu))^ and k^n"[ lu) is 
the attenuation factor for a surface plasma wave (wherein 
electron-level damping is included). The continuous field 
operators obey [c(a;, x), , x')] = 5{x — x')5{uj — oj') 

and the second term in Eq. \27\ preserves the bosonic 
nature of the propagated SPP. 

Using the relation b{t) = (27r)“^/^/ da;e“*‘^‘ 6 (a;), the 
mean SPP flux at space-time coordinate (x, t) can be cal¬ 
culated, foutix,t) = {b\^^{x,t)bout{x,t)). Eor a narrow 
wavepacket centered at wq, we have [33] 


/out(x,t)=e- 2 '=o«"Vout(ifl), (28) 


where Ir = t — x/vg, with Vg being the group velocity 
at the center frequency. The mean flux of the quantized 
SPPs is therefore simply damped by the classically ex¬ 
pected factor 2kQK". Using the values of k" in Pig. for 
both TE and TM modes, and the above model for the op¬ 
erator mapping (summarized by Eq. (Hii)) we are now in 
a position to further investigate the impact of loss on the 
transfer of quantum states of SPPs on the graphene sur¬ 
face and quantify the performance of an error-correction 
code for protecting against this loss. 


V. ROBUST-TO-LOSS QUANTUM STATE 
TRANSFER 

A. Lossy propagation 

In Section III we showed that efficient coupling of inci¬ 
dent single photons and graphene SPPs is possible at the 
quantum level. In this section we now analyze the trans¬ 
fer of more complex photon states to SPP states, and 
their subsequent propagation. The general setting is the 
following, the input state is defined as |4')in = \'tp}a\0)b, 
where a corresponds to the photon mode and b to the 
SPP mode (which is initially in the vacuum state). The 
photon interaction with the SPP via end-fire or prism 
coupling produces the output state, written as |'I')out = 


77|^')in, where the unitary transformation U is defined by 
Eq. ( |^ . We start by considering the photon input in a 
superposition of coherent states [33], 


|^')i„ = A(|a) + |-a)),|0)b, 


(29) 


with I ±a)a = exp[-|ap/ 2 ] where 

|n) is a number state and N = {2 + Us¬ 

ing the transformation matrix from Eq. (26 1 we have 
|4')out = A^(|acos 5 )a| -asing)b + \ - a cos^|a sin^)^) 
where j5 = sing and 7 = cosg. For perfect coupling 
{g = 7r/2) we have 


|«')out = |0)aiV((|-a) + |a))b (30) 


and the coherent state superposition is transferred per¬ 
fectly to an SPP superposition. 

In general we are interested in the SPP state itself, 
so we need to trace out the unobserved photon modes 
a. The density operator for the total photon and SPP 
system is p = |'I')out out(^'|- By tracing out system a we 
have pb = Tr^p, which gives [33] 

Pb = |fVp(|asing)(asing| -|- | — asing)(—asingj 

-I- Co„t(|asing)(-asing| -|- | - asing)(asing|)) (31) 


where Cout = exp[—2|acosgp]. We take this state as the 
initial mixture (at x = 0 ) that we want to propagate a 
distance x along the graphene in the presence of loss and 
characterize its decoherence. Using Eq. (271 and Eq. (31) 
one finds mill] 


P^{x) =|iVp (1 - asinge-'=«”"^)(-asinge-'=“'^"^| (32) 
-f |asinge-'=“«'^)(asinge-'=»'="“| 

-h c(x)(|a sin ge-^^^"^){-a sin ge-^'>^''^\ 
+ I - asinge-*^“'‘"“)(asinge-'=«'""“|)) 


where c(x) = Cout exp[—2|a singp(l — “)]. Note 

that at long times (large x) the SPP moves towards the 
vacuum state, as expected, and at early times (small x), 
c(x) ~ Cout, and therefore Pb(x) ~ Pb(0). 

In Section III we showed that good coupling can be 
achieved using several different coupling methods (end- 
fire, which has been experimentally demonstrated [44], 
and prism coupling [35H1I]), resulting in values of g ~ 0.8 
or higher. In the following investigation of state propa¬ 
gation, rather than link the results to a certain coupling 
geometry we assume values of g in a reasonable range. 

In order to quantify the effect of loss on the excita¬ 
tion process and subsequent propagation, we use the fi¬ 
delity F = {'\j)\p\'ip) as a measure of the similarity be¬ 
tween two states, one pure (ideal) and one mixed p 
(damped) [52] . When F = 1 the states are the same and 
when F = 0 they are completely orthogonal. In the ideal 
case, the superposition state \tj}) given in Eq. (29) will be 
excited as an SPP and then propagate without loss along 









FIG. 5: Fidelities of propagating SPPs with respect to the 
initial state in (a) and the initial state with smaller ampli¬ 
tude OL = asinge“*^“” ^ in (b) after excitation with differ¬ 
ent excitation couplings /3 = 1,0.98,0.95,0.9,0.8. The initial 
photon state has 0 = 3. The plots in (a) monotonically de¬ 
crease to zero, while those in (b) approach a finite value when 
koK"x < ln[%/2/3], after which they turn and increase to unity 
as the vacuum contribution becomes significant. 


the graphene surface. In the realistic case, however, the 
state p given in Eq. (321 will be the state resulting from 
the non-ideal excitation process and damping. The fi¬ 
delity thus provides a means to measure how far away in 
the Hilbert space the damped SPP state is from the ideal 
(initial photon) state as it propagates along the graphene. 

In Fig.[^(a) we show examples of the fidelities between 
the initial photon state in Eq. (291 and the excited SPP 
states in Eq. 


(with different coupling efficiencies) 
which become damped as they propagate. One can see 
that the initial fidelities do not start at 1 for the non¬ 
ideal coupling cases and the fidelities decay as the SPP 
propagates, showing the movement of the quantum state 
further away from the initial ideal state. In Fig. [ 11 (b) 
we show the fidelity of the excited SPP state with re¬ 


spect to the ideal photon state with a smaller amplitude 
a' = as\Tvge~^°'^ We have included this case as the 
amplitude of the SPP state is expected to decay as it 
propagates, thus it is informative to compare it with an 
ideal photon state that has a decayed amplitude, but 
importantly has no degradation in its original structure 
- it is a pure state with the same fixed positive phase. 
In this second case the fidelity starts at a slightly higher 
level and approaches a higher asymptotic value (= 0.5) as 
the damping increases. In Fig.[^(a) the asymptotic limit 
of the fidelity is zero as the SPP state moves toward the 
vacuum state as a result of dissipation of energy. On the 
other hand, in Fig.[^(b) the asymptotic limit is 0.5 as the 
state we are comparing the SPP state with is matched in 
terms of its energy, resulting in an effective phase damp¬ 
ing of the SPP state from the perspective of the photon 


state. 

From Fig. it is clear that the excitation process and 
subsequent damped propagation affect the quality of the 
quantum state transfer between photons and graphene 
SPPs. In the next section we will show that by using an 
error correction strategy, one can protect the superposi¬ 
tion state (and more general quantum states) from loss 
caused by the damping during propagation. 


B. Error-correction code 

In order to provide robust-to-loss propagation of quan¬ 
tum states of SPPs along the graphene waveguide we 
consider the following code states [55| 


| 0 ±) = 

^(H±|-a)), 

(33) 

|i±) = 

y^(l*a)±l *«))> 

(34) 


where N± = 2(1 ± ). For large enough mean pho¬ 

ton number (h) = |ap we have that the states |±a) 
and \±ia) are orthogonal, and therefore so are the code 
states. The states |0+) and |T+) form an orthogonal ba¬ 
sis, |(l+| 0 +)p « 0 when a > 2 , representing a code 
space, in which an arbitrary quantum bit (qubit) can be 
encoded as 


Id/) =co|0+)+ci|T+), (35) 

where |coP + |cip = 1. This type of encoding has re¬ 
cently been considered theoretically in a cavity scenario 
using superconducting qubits [siins] and experimentally 
demonstrated in Ref. [56]. Here, we extend its applica¬ 
tion to a waveguide setting. 

When damping occurs during SPP propagation, for 
a short time interval At this can be described by the 
density matrix mapping m 

P ^ P — Pjump “1“ (1 Pno—jump-} (36) 

where AP ^ I is the probability of losing an 
excitation during time At, AP = 7 AtTr(/ 3 h) = 
jAt{h), with 7 the decay rate. Here we have 
more explicitly Pjump = apa^ {fi) and Pno-jump = 

g-JAta+a^gJAtata/rp^Jg-JAiata^gJAta+aj^ 

words, we can write the evolution of the density matrix 
during a short time At as p —>■ p' = EqpE\ + EipE\, 
where the Kraus operators for the damping process are 
Eq = and Pi = (I - ^Ata^'a), with ~ 

(1 — ^Ata^'a). By using the relations 

a|0±) =ax/iV^/iV±|0^), a\l±) = iay'N^/N±\l^){37) 
g-3Atata|J^^ ^ gll“l^(<=-"^‘-l),yiV±,At/fV±|I±.At), 
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where iV±_At = 2(1 ± e ^ ^ I ) and the decayed 
states 

|0±,At) = (38) 

|I±,At) = ^^(|*ae-'^^‘/^)±|-*ae-^^‘/2)),(39) 

y/^^±.At 

one finds that the code states are mapped as follows 

|0+)(0+| ^ AP|0_)(0_| + (l-AP)|0+,At)(0+,At| 

|I + )(I+| ^ AP|I_)(I_| + (l-AP)|I + ,At)(I + ,At|, 

and similarly for the off-diagonal terms. The code space 
and the decayed code space {|0+,At), |i-i-,At)} - which re¬ 
mains orthogonal for large enough mean excitation num¬ 
ber, can be distinguished from the erred space (spanned 
by |0_) and |1-), and their decayed versions) by the 
photon parity operator P = e“a’^a = = 

l)”|n)(n|. Using the relations Pla) = |—a) and 
P|—a) = |a), one finds (0+|P|0+) = (i+|P|l+) = +1 
and (0_|P|0_) = (I_|P|I_) = —1. Thus, by measuring 
the parity continuously within small enough time periods 
one can determine whether or not an excitation has been 
lost and correct the state back into the code space. 

Two important points should be mentioned in relation 
to our application of the above error correction strategy 
to propagating SPPs in a waveguide setting. First, the 
correction operations do not need to be performed until 
the very end of the propagation. This is because the par¬ 
ity checks continuously project the state into either the 
code space or the erred space, with the state moving be¬ 
tween these two subspaces in a cyclical fashion as it prop¬ 
agates, similar to the case described in Refs. [MUSS]. At 
the end of the propagation, after having recorded the out¬ 
comes of the sequence of parity checks, we know whether 
the final state is in the code space or the erred space and 
we can correct it accordingly to bring it back into the 
code space, although with a reduced amplitude. This 
leads to the second important point, which is that once 
the initial state has propagated a given distance (and un¬ 
dergone many parity check operations), it will have de¬ 
cayed significantly. Therefore the error correction strat¬ 
egy means that we must supply a large enough starting 
value of a for a desired propagation distance along the 
graphene surface so that the decayed code states main¬ 
tain their orthogonality. By doing this we effectively 
put quantum state transfer and classical state transfer 
in plasmonics on a level playing field, where one simply 
increases the intensity in order to transfer information 
along the plasmonic waveguide. 

In summary, the robust-to-loss encoding for quantum 
state transfer requires only parity checks to be performed 
on the SPP as it propagates along the graphene sur¬ 
face. These checks could be realised by incorporating 
additional circuitry within the graphene sheet, using an 
ancilla mode (electronic or photonic in origin) and the 
unitary operation outlined in Ref. [56j . which provides 
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FIG. 6: Average fidelities of error-corrected propagating 
SPP qubits with increasing parity-check probabilities p = 
0,0.2, ■ • ■ , 1 for different excitation couplings: /3 = sin[^ = 
7r/2] in (a), /3 = sin[(ji = 0.9 x 7r/2] in (b), and /3 = sin[g = 
0.8 X 7r/2] in (c). Fq is the average fidelity just after the ex¬ 
citation process, before the first parity check operation. The 
fidelities are calculated with respect to the initial state with 
a smaller amplitude a' = as'mge~'‘°'^ ^ for a = 3. The prop¬ 
agation length is scaled with respect to TM (TE) modes at 
T=300 (0) K and pc = 1.4 (0.4) eV for Ao = 1550 nm in 
top (bottom)-upper horizontal axis and Ao =810 nm in top 
(bottom)-lower horizontal axis. The effective wavelengths for 
TM (TE) modes for Ao = 1550 nm and Ao = 810 nm are 
Aeff(= 2TvlkoK') ~ 36.42 (1549.77) nm and ~ 7.23 (810.053) 
nm, respectively. Panel (d) checks the orthogonality approx¬ 
imation, |(0±|1±)I < 10“^, is well satisfied during a quantum 
jump simulation for the case of (c), for instance. 


a means of carrying out a non-invasive measurement of 
the photon parity. The closeness of the parity checks 
on the surface depends on the rate of loss of the SPP 
and the speed with which it propagates. For this, we as¬ 
sume the SPP wavepacket containing the quantum state 
is narrowly centred around the frequency Wq, such that 
it propagates with speed vci^Jo) = vq- We then use 



— I<1.'|0/>l 

— I<1-'|0_'>| 
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the relations 7 = 2koK.''va and At = Ax/vq to find the 
corresponding distance Ax for the time interval At. 

In the following analysis of the performance of the 
code, we carry out a quantum jump simulation via a 
Monte Carlo iteration method, tracking 10^ individual 
trajectories of the initial state and summing the final 
states m- We use the fidelity to quantify the effec¬ 
tiveness of the code as the number of parity checks 
is modified from the ideal case when it occurs every 
Ax <C {2koK''\a\‘^) m, to the case where no parity checks 
are made, corresponding to bare graphene propagation. 
In all cases we consider an initial parity check after the 
photonic excitation of the SPP in order to put the differ¬ 
ent cases on an equal footing from the point of excitation. 
As examples, we show the fidelities (averaged over the 
single-qubit Bloch sphere) of error-corrected propagation 
of TM (TE) modes at T=300 (0) K and = 1-4 (0.8) 
eV for Aq = 1550 nm and Aq = 810 nm in Fig. for 
different excitation couplings in (a), (b), and (c). During 
propagation, the parity check is performed every At/p 
(equivalently Ax/p) for a given parity-check probability 
p. Note that the proposed scheme corrects the damp¬ 
ing effect by flipping over the states in the erred space 
depending on the outcomes of the parity measurements, 
leading to a significant increase of fidelity from the initial 
value Fq given by the excitation process. In Fig. §(d) 
we show the validity of the orthogonality approximation 
of the code basis states for the fidelity calculation as the 
SPP propagates, which we set via |(0(j_|l±)| < 10“^. 


VI. CONCLUSIONS 

In this work we investigated the excitation efficiency 
and impact of loss on SPP propagation in graphene at 
the quantum level. We considered two different excita¬ 
tion techniques: end-fire and prism coupling. We started 
by quantizing the transverse-electric and transverse- 
magnetic SPP modes in graphene, and used this to build 
a fully quantum model for the excitation process. We 
then studied various parameter regimes that enabled the 
excitation of SPPs by photons and found that efficient 
coupling of single photons to graphene SPPs is possi¬ 
ble. We then studied the subsequent propagation of ex¬ 
cited quantum states under the effects of loss induced 
from the electronic degrees of freedom in the graphene. 
In order to protect the quantum states from loss we 
used a quantum error-correction code and found that the 
code provides a robust-to-loss mechanism for propagat¬ 
ing quantum states of light in graphene over large dis¬ 
tances. The results and analysis in this work contribute 
to the growing field of quantum plasmonics, and to the 
use of graphene as a flexible alternative to basic metallic 
materials for supporting SPPs and their quantum appli¬ 
cations. 


0.76 



0 . 62 '-^^^^^- 

9 9.5 10 10.5 11 11.5 12 

Frequency (THz) 


FIG. 7: Grating end-fire coupling of photons to SPPs on 
graphene. Inset shows the configuration considered and the 
main part shows the transmission coefficient p. The grating 
is silicon, has period 4.717 ^m and depth 2.3 pm. The width 
(across) of the graphene region is 39.5 pm, and pa = 0.152 
eV, r = 0.1 ps at r = 0 K. 
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APPENDIX 

Appendix I: Pree-Space Photon to SPP Coupling 
Models 

In the main text we considered photon-to-SPP quan¬ 
tum state transfer for several values of the coupling pa¬ 
rameter g. In this appendix we provide several coupling 
models that can be used to achieve good coupling. 


A. End-Fire Coupling 

In the first configuration, we consider end-fire coupling 
from a fiber (TE) or wire (TM) to graphene strips. For 
the TE case it has been experimentally verified that good 
fiber-to-graphene coupling exists |44j . THz experiments 
for the TM case have not been performed as yet, although 
far-infrared excitations of TM SPPs has been shown us¬ 
ing an atomic force microscope tip [SSj [SS] . 

For the TE case we consider a structure similar to the 
experimental configuration of Ref. [44) . which consisted 
of an optical fiber (step-index core-cladding) with a cen¬ 
ter perturbed region where the cladding is removed and 
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the core partially removed, upon which the graphene is 
located (see inset to Fig. |^(a)). The fiber core is excited 
at the far-right as the input and the output is taken at 
the plane on the left where the structure terminates. In 
Fig. we show the transmission coefficient j3 obtained us¬ 
ing FDTD simulation via CST Microwave Studio. Very 
good selectivity for TE polarization is exhibited, as mea¬ 
sured in Ref. [H]. Note, the dimensions of the structure 
in Ref. [U] and those used here are somewhat different 
(see caption for details). Good TE SPP propagation is 
expected in the frequency range shown, since here a" < 0 
and only a TE SPP propagates. Fig. [^(b) and Fig. [^(c) 
show the power distribution on the structure for TE and 
TM excitation respectively, showing good TE SPP prop¬ 
agation and no TM SPP propagation. 

For the TM case at low THz frequencies, a grating 
geometry can be used for SPP excitation. The inset of 
Fig.0 shows the coupling geometry and the main part 
shows the transmission coefficient /3. At these frequencies 
a" > 0 and only a TM SPP can propagate. 

For both the fiber and grating end-fire coupling meth¬ 
ods, the quantization of the graphene SPPs follows the 
same method as that given in the main text, with appro¬ 
priate consideration of the asymmetry in the surrounding 
dielectric media. 



1 1.5 2 1 15 2 

Freq. (THz) Freq. (THz) 


FIG. 8: Reflectance for the prism-coupling configuration, (a): 
Reflectance i? = as a function of frequency and incident 
angle for the TE case, for a prism having £i = 1.5, prism- 
graphene spacing d = 620 nm, and /Tc = 1.24 eV. (b): TM 
reflectance for si = 1.5, d = 200 pm and pc = 0.5 eV. (c): 
TM reflectance for ei = 3.9, d — 200 pm and pc = 0.5 eV. 


B. Prism coupling 

For prism coupling of photons to SPPs we consider 
an attenuated total reflection (ATR) set-up (Otto con¬ 
figuration |60| 1 to provide a momentum match between 
the incoming photon and the graphene SPP. We ignore 
photon reflection from the top of the prism (which can be 
mitigated using impedance matching), and assume that a 
plane-wave photon field is incident on the lower air-prism 
interface at z = d, as depicted in Fig.|^(a). Total internal 
reflection will result in an evanescent field in the space 
below the prism, which can couple to the SPP evanescent 
field. Since the space below the graphene is vacuum, the 
field transmitted into that region is also evanescent. For 
graphene, prism coupling to SPPs has been considered 
classically in Refs. [45H47j. 

To determine the efficiency of SPP excitation we con¬ 
sider the overlap between the photon field and the SPP 
field. For the photon field we assume an incident field 
having parallel polarization. In Appendix II we quan¬ 
tize the TE and TM photon modes for the ATR geom¬ 
etry, leading to the quantized vector potential operator 
for photons having the form 

Ap(r,t)cx J dw(r(a;)T'“(r,a;) 0 ( 2 ; — d) 

-I- T{u!)'t>^{r,uj)Q{d - z))d{uj)N^^^‘^{uj) + H.c., 

where Nqr is a normalization parameter having units of 
length, and r{uj) and t{uj) are field amplitudes. The fac¬ 
tors are mode functions that depend on the 


geometry. The upper-region mode function T'" has a 
standing wave behavior in z and cannot couple to the 
graphene plasmons, so that this term can be ignored. The 
lower-region mode function T'^(r,a;) = T'™(r,a;)(0(z) — 
0(z — d)) + 'l'‘(r,a;)0(—z) exhibits evanescent behavior 
and couples energy into the SPP. 

As described in Refs. [311113], the transmission coeffi¬ 
cient is the overlap integral between the two fields (pho¬ 
ton and SPP), given by 


P*{u}) = a/ 1 - |r(a;)|2(5(a; - u}')d{k - kl) 


(40) 




= y dzT'(z,w) • </.(z,w')*. 


Here, uj {oj') is the radian frequency of the photon (SPP), 
kl. (fc) is the X component of the photon (SPP), is the 
SPP normalization and the SPP mode functions 4>{z, uj') 
are given explicitly in Eq. The form of the ATR 

mode functions is given in Appendix II. The transmis¬ 
sion coefficient j3 depends on the geometrical parameters 
of the prism-graphene system (governing the reflection 
coefficient r and the degree of mode overlap), and repre¬ 
sentative results are presented in the following discussion. 


In order to excite SPPs on the graphene surface, the 
longitudinal wavenumber of the incident field, kl = 
kl sin^i, must match the propagation wavenumber of the 
SPP given in Eqs. (|^ and (|9), where ki = y/elko and Si 
is the relative permittivity m the prism. Therefore, set¬ 
ting = kl sin 9i leads to the matching frequency ojq 
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that satisfies 


1. TE prism modes 


a{oJo)'r] = ±2^ 1 — £i sin^ 9 


a{uJo)ri = 


±2 




£i Sin 


TE, 

TM, 


where 0 > 9^, the critical angle for total internal reflec¬ 
tion (TIR). However, this 9 is only real-valued for lossless 
graphene (ct = ia”). For the more realistic lossy case one 
must find a zero or minimum of the reflection coefficient 
r at the prism-air interface. 

In Fig. (a) we show the reflectance R = |rp as a 
function of frequency and incidence angle for the TE case, 
for a prism having Si = 1.5, prism-graphene spacing d = 
620 nm, and /ic = 1.24 eV. The critical angle for TIR 
is 9c = sin“^(l/Y^) = 54.736° (for 9 < 9c the SPP is 
not excited). Since for the TE mode, which is loosely 
confined to the graphene surface, ~ fco, a match is 
found for ei sin^ 0 ~ 1, which occurs very close to the 
critical angle. As a result, the contour of the coupling 
angle is a horizontal line near 0 ~ 0c- 

For the TM case at low THz frequencies. Fig. |8](b) 
shows the reflectance for £i = 1.5, d = 200 /rm and 
/ic = 0.5 eV, and Fig. (c) shows the result for a prism 
having Si = 4. Compared to the TE case, considerable 
dispersion is found, with good matching at low THz fre¬ 
quencies and angles moderately above the critical angle. 

In Fig. |^(b) we show the transmission coefficient /3 as 
the prism-graphene spacing d changes for TE (upper hor¬ 
izontal axis) and TM (lower horizontal axis) SPPs. The 
inset shows the TE reflection coefficient versus frequency 
near resonance, which occurs for 9 9c and fku ~ 2^c- 
Although the inset only shows the TE case, the reflection 
coefficient can be reduced to zero for both the TE and 
TM cases. Despite this, the overlap integral results in 
/3 < 1 since the mode functions outside the prism region 
are not identical to those inside the prism. Since we do 
not account for reflection from the prism-graphene inter¬ 
face in a rigorous manner, our calculation of beta is valid 
when the EM field energy inside the prism is negligible. 


For the TE case (perpendicular polarization), we as¬ 
sume an incident field in the prism (z > d) given by 

E*(r) = -h c.c., 

kl-k - kiz 


H*(r) = 


UfJ-i 


. + C.C., 


where fcj, = fcisin0i, fc® = fcicos0i and ky = 0. The 
reflected field in the prism is 

E’’(r) = -h c.c. 

^ klx - klz i{klx+k:z)^-iut 

with fcO = fc® and kl = — fc®. The field in the middle 
region (below the prism and above the graphene, 0 < 
z < d) is 

E'"(r) = -I- 

-I- C.C., 

l,ml~ _ J,ml* 

r\j ^ A. ^ 


H’"(r) = 


w/ri 
r,m2y 

IVfii 


^ j^ml ^i{k^x-£k^ z) 


£;m2gqfc™x-/crDe-“‘ + C.C., 

(41) 

where A™ = k^. and fc™ = — (/c^,)^. The transmitted 


field {z < 0) is 


E*(r) 

H*(r) 


-Pc.c., 

~ EW{kix+klz)^-^UJt 


with fc* = and fc* = fc™. Enforcing the boundary con¬ 
ditions of Eq. Q, we obtain the reflected and transmit¬ 
ted field amplitudes E]_, E^, E™ and The resulting 
vector potential has the form 


A(r,a;) 

= r(w4'“(r, t)Q{z — d) + r(a;)4''^(r, t)Q{d — z) + c.c., 


where 4'“>^(z,a;) are mode functions, and r{uj) = E^/E\_ 
and t{uj) = E'\_ are reflection and transmission 

parameters. 


Appendix II: ATR Fields and Quantization 


2. TM prism modes 


To compute the reflection coefficient in Eq. (40) for 


the ATR geometry, we assume an incident photon field 
and solve the plane wave reflection/transmission problem 
for the prism-graphene geometry shown in Fig. We 
then quantize the resulting fields and obtain the mode 
functions. 


For the TM case (parallel polarization), we assume an 
incident field 

E'^ 

H*(r,t) = -hc.c.. 


E*(r,t) = 


m 

-klx + klz 


kl 


c.c. 
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where = fcisin^i, fc® = ki cos 9i. The reflected field 
{z > d) is 


E'^ 

H'’(r,t) = + C.C., 

»?i 

+C.C., 

fci 


where are mode functions, and r(a;) = i?y/£iy 

and t{u}) = A/fc‘/fc*e®^*‘^i?||/i?y are reflection and trans¬ 
mission parameters. 


3. TE and TM prism mode quantization 


with k^. = kl and k^ = —k\. The field in the middle 
region (below the prism and above the graphene, 0 < 
z < d) is 


For either the TE or TM polarization the Hamiltonian 
can be evaluated in the same way as for the graphene 
quantization described in the main text, leading to 


H™(r,t) 

E™(r,f) 


E^ ^m2 

(_Lg-i/cr^ + -h C.C., 

m Vo 

_U'rn- I l.m- 

i%z ^ ^ml ^i{k^ x+k^ z) 


Um- I h,m- 


H = eoL^ Y. + <^L), (42) 

k 

where iVk has units of length, so that the field is quan¬ 
tized using 


with fc™ = and fc™ = yjk^ — {kl)"^. The transmitted 
field {z < 0) is 


H‘(r,t) 

E‘(r,t) 


E^ 

+C.C., 

Vo 

— fc(,X-|-fc*Z i{ktx+kiz) -iujt 

ko II 


with fc* = and /c‘ = fc™. Enforcing the boundary 
conditions of Eq. 0 , we obtain the reflected and trans¬ 
mitted field amplitudes. The resulting vector potential 
has the form 


A(r,a;) 

= (r(a;)4'”(r, t)0(z — d) + T{uj)'$^{r, t)Q{d — z)) + c.c.. 




At^ 


h 


2soL^ujNi, 

h 


^k 5 


(43) 


(44) 


2eoL‘^u}N^ 

The bosonic annihilation and creation operators dk and 
dk satisfy [dk,dk,] = dk.k'- Assuming a continuum of 
frequencies, the quantized potential operator is then 


Ap(r,t) = 


1 


H 


(45) 


2tt J ■■ V 2eoL'^ujN 

X (r(w)4'“(r,w)0(z-d) 

-I- r(tLi)'Ic^(r, ui)Q{d — z)) + H.c. 
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